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SECTION  I. 


INTRODUCTION 

1.  Objective  of  the  Study 

In  this  report,  analytical  models  are  formulated  for  obtaining  the 
probability  of  kill,  P^,  of  a  point  target  for  a  variety  of  situations.  These 
situations  can  be  roughly  divided  into  two  categories,  namely,  multiple 
weapons  delivery  (Sections  II  and  III)  and  single  guided  weapon  delivery 
(Section  IV  and  V). 

In  Section  II,  the  general  problem  of  the  stick  delivery  of  weapons  's 
formulated.  An  exact  expression  for  is  obtained.  However,  this  expression 
is  in  the  form  of  a  multiple  integral  which  is  very  dificult  to  evaluate.  As 
a  consequence,  an  approximation  is  suggested  and  the  reader  is  referred  to 
for  further  details. 

In  Section  III,  the  problem  of  multiple  weapons  delivery  with  independent 
passes  is  described.  In  general,  at  each  pass,  weapons  are  delivered  in 
stick,  and  thus,  a  formal  approach  to  the  problem  would  have  to  rely  on  the 
results  of  Section  II.  No  specific  methodology  is  provided  to  tackle  the1 
general  problem.  The  only  situation  that  is  worked  out  in  detail  is  the  case 
of  singly  delivered  weapon  per  pass. 

Sections  IV  and  V  consider  the  problem  Of  guided  weapon  delivery.  The 
delivery  of  guided  weapons  differ  from  the  delivery  of  general  purpose  bombs 
In  that  the  delivery  error  distributions  have  different  characteristics. 

In  its  simplest  form,  for  guided  weapons,  the  delivery  error  in  the 
direction  of  range  has  the  same  normal  distribution  as  for  the  delivery  error 
in  the  direction  of  deflection.  The  two-standard  deviations  are  equal.  Given 
that  the  two  delivery  errors  are  unbiased  and  that  they  are  independently 


distributed,  it  immediately  follows  that  the  distribution  of  errors  in  the 
radial  direction  follows  the  so-called  circular  or  Rayleigh  distribution. 
Experimental  evidence,  however,  has  shown  otherwise  and  two  attempts  are 
made  In  the  JTCG/ME  (Alr-to-Surface)  manuals  to  correct  for  this  situation. 

In  the  early  edition  (1980)  [1],  a  modified  radial  distribution  is  proposed 
for  laser-guided  bombs.  This  situation  is  studied  In  detail  in  Section  IV. 

In  a  more  recent  edition  (1983)  [2],  a  normal  delivery  error  with  bias  In  the 
range  direction  is  proposed  for  general  guided  weapons.  This  situation  is 
studied  in  Section  V.  Section  VI  consists  of  concluding  remarks. 

2.  Assumptions 

The  following  -ssumptions  are  made: 

a.  Each  of  the  target  and  weapon  is  idealized  as  a  point.  Depending  on 
the  particular  situation  considered,  the  weapon  may  or  may  not  be 
aimed  directly  at  the  target. 

b.  The  direction  of  the  weapon  delivery  range  and  deflection  are 
respectively  parallel  to  the  (x-y)  coordinate  system  on  the  ground 
plane.  Since  the  coordinate  system  can  be  arbitrarily  selected,  there 
is  no  loss  in  generality  in  making  this  specific  assumption.  The 
position  of  the  target  has  coordinates  (u,v). 

c.  For  stick  delivery,  the  weapons  are  subject  to  ballistic  and  aiming 
errors  in  each  of  the  x  and  y  directions.  For  single  weapon  release, 
a  delivery  error  in  each  of  the  x  and  y  directions  is  assumed.  In 
each  case,  the  probability  density  function  of  these  errors  wist  be 
specified. 


Z 


The  probability  of  kill  due  to  fragmentation  at  a  point  (u,v) 

given  that  the  weapon  Impacts  at  (x,y)  is  given  by  the  three-parameter 

Carleton  damage  function 

D(u-x,  v-y)  *  l)0  exp{-  D0[tnr)2  + 

x  y 

'ne  pnrameters  Rx  and  Rv  are,  respectively, the  weapon  radii  in  the  x 
and  y  directions.  The  parameters  Dq  is  the  maximum  probability  of 
kill  which  occurs  at  the  point  of  weapon  impact 
Fragmentation  does  not  contribute  to  the  delivery  error. 
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SECTION  II 

MULTIPLE  WEAPC,  iS:  STICK  DELIVERY 


1.  Introduction 


We  consider  the  stick  delivery  of  n  weapons  (n»l,2,...)  for  which  the 
following  assumptions  are  made: 

a.  Each  weapon  is  subject  to  ballistic  errors  which  are  assumed  to  be 
normally  distributed,  independent  of  each  other  and  independent  in 
each  of  the  x  (range)  and  y  (deflection)  directions. 

b.  The  entire  stick  pattern  is  subject  to  aiming  error  which  is  assumed 
to  be  normally  and  independently  distributed  in  each  of  the  x  and  y 
directions,  and  independently  distributed  from  the  ballistic  errors. 
Equivalently,  it  can  be  stated  that  each  of  the  n  weapons  is  subject 
to  the  same  aiming  error. 
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Figure  1.  Geometry  of  Target  and  Weapon 


Let  (C1 >  n^)  be  the  mean  point  of  Impact  (MPI)  of  the  ith  weapon 
(1*1,2,... ,n).  The  ballistic  errors  in  the  directions  of  range  and  deflection 
are,  respectively,  and  Yj.  The  random  variables  X1  and  Yj  are  assumed  to 
be  Independent  each  with  zero  mean  and  having  the  respective  normal 
probability  density  functions 


For  each  of  the  i  weapons,  the  aiming  errors  in  the  directions  of  range 
and  deflection  are,  respectively,  X  and  Y.  The  random  variables  X  and  Y  are 
assumed  to  be  Independent  each  with  zero  mean  and  having  the  respective  normal 
probability  density  functions 


9 v(x)  •-TTi:  exp[--*r]  (3) 

a/2.  2ax2 

and 


.  2 

gY(y)  * — exp{--^]  (4) 

•/ST  to* 

The  random  variables  { } ,  X,  and  Y  are  actually  Independent. 

Let  (W^ ,  Z^)  be  the  actual  Impact  point  of  the  1th  weapon.  Thi3  is  the 
result  of  the  combined  effect  of  the  ballistic  error  and  aiming  error.  This 
combined  effect  is  the  sum  of  the  ballistic  error  and  aiming  error  in  each  of 
the  x  and  y  directions  (see  Figure  1).  Clearly,  fpr  1«l,2,...,n 


“1  «  (S) 

*1  ■ ",  *  V  ’  !S) 

The  sequence  of  random  variables  {W^}  and  {Zj}  are  "wtually  Independent  and 
are  normally  distributed  with  respective  means 

ECW-1  *5^  ECZ.'J  -  n1  (7) 


and  respective  veriances 

Var[W.]  =  +  a*;  Var[Z.1  *  o*  +  °y  v*) 

The  random  variables  in  the  sequence  {W,.},  1*1 . .  are  not  mutually 

Independent  since  all  of  them  depend  on  the  common  random  variable  X.  Clearly 
for  i*j 


Cov[W. ,  W-]  =  EKw^Xw.-^)] 

*  E[(Xi+X)(X*+X)l 

=  E[X2]  »  oZx  (9) 

Thus,  the  sequence  of  random  variables  {W^|  are  jointly  normally  distributed 
with  mean  vector 

*  (c  ^  2  *  *■*  • 
and  variance-covariance  matrix 


cz 


6 


Thus,  the  sequence  of  random  variables  {Z ^ ,  1*1,2,... ,n  are  jointly  normally 
distributed  with  mean  vector 


n_  a  (Oj,  »nn) 


(15] 


and  variance-covariance  matrix 


7  7 

°2  +  ay 


A  = 


2  .  7 

a2  +ay 


I  n2  2 

I  •  O  0  • 

L  y  .  y 

2  2 
=  I  *  ri 

If  we  define  Z_  *  (Z1,Z2,...,Zn) ,  then 
of  the  Z | * s  is  N(Z_;  n,  a). 


2  2 
a2  +  ay 


the  joint  probability  density  function 


2.  The  Model 

For  a  target  located  at  (u.v).  It 
bility  of  kill  of  the  target  assuming 
damage  function  Is  of  the  form 


(16) 


(17) 


Is  required  to  determine  the  proba- 
hat  If  a  weapon  Impacts  at  (w ,z),  the 
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(18) 


0(u-w,  v-z)  =  D0  exp{-  Og^-g-)  +  HH  M 

x  y 

The  conditional  probability  of  kill  given  that  the  first  weapon  impacts 

at  (wj.Zj),  the  second  weapon  impacts  at  . and  the  nth 

impacts  at  (wn,zn)  is 

{1  -  n  [1  -  D(u-w. ,  v-z,)]}  (19) 

i*l  11 

The  unconditional  probability  of  kill  at  (u,v)  is 

OS  00 

P.  =  /  ...  /  [Probability  of  kill  at  u,v)  |  ith  weapon  impacts  at  (w^,z^) 

2n 

for  i>l,2,...,nl[Probability  that  the  ith  weapon  impacts 
between  (w^.Zj)  and  (wj+dwj,  z^+dz^)  for  i=l,2,...,nl 

m  m  n 

*(.../  { 1  -  H  [  1  -  D(u-w^ ,  v-z^]}  N(w;  £,  ft) 

-Tr' 

n 

N(z;  n.  A)  H  dw,  dz^  (20) 

-  1*1  1  1 

Note  that  P^,  in  general,  is  going  to  be  a  function  of  the  location  of 
the  target,  namely  (u,v),  and  also  a  function  of  the  cordinates  of  the  MPI  of 
all  n  weapons.  It  is  customary  to  assume  that  the  location  of  the  taget 
coincides  with  the  origin,  so  that  u«0«v,  and  hence,  P^  is  only  a  function 
of  £  and  n_  and  one  can  write: 

P.{£.n)  =/.../  {1  .  I  [1  -  0(w1  N(w,  L,  ft) 

*  i*l 

2n 

n  • 

N(z,  n.  A)  a  dw,  dz,  ^l] 


3.  Remarks 


There  are  two  problems  to  be  addressed  here.  The  first  problem  is  that 
of  obtaining  the  optimum  stick  pattern  that  will  maximize  the  probability  of 
kill  P|<U,n).  By  an  optimum  stick  pattern  is  meant  the  determination  of  the 
optimum  locations  of  the  MPI's  of  each  of  the  n  weapons  delivered  in  such  a 
way  that  the  probability  of  kill  is  maximized.  As  one  of  the  possible  courses 
of  action  one  could  select  the  MPI's  to  coincide  with  the  location  of  the 
point  target,  that  is  to  set 

£ j  *  0  *  ,  1  *  1,2,... ,n 

But  this  course  of  action  does  not  necessarily  maximize  the  probability  of 
kill.  Note,  however,  that  the  weapons  are  dropped  from  specific  locations  on 
the  plane.  There  are  usually  three  such  locations:  the  two  wings  and  the 
airplane  centerline.  As  many  as  100  of  these  weapons  may  be  released  In  a 
single  stick  delivery.  Thus,  from  a  practical  point  of  vlrw,  the  variables 
and  are  selected  according  to  a  constrained  path. 

The  second  problem  consists  in  obtaining  a  computable  expression  for 
P^i.  H)*  2n-  tuple  Integral  could  theoretically  he  evaluated  (since  the 

function  n(w^,  z^)  is  of  an  exponential  form).  However,  It  is  extremely 
difficult  to  proceed  with  the  Integrations  particularly  for  large  number  of 
weapons.  The  complexity  of  the  expression  for  has  led  to  the  develop¬ 

ment  of  an  alternate  method  for  rederiving  PjJi.nj.  The  other  method  relies 
on  a  decomposition  principle  v^lch  results  in  a  substantially  more  simplified 
expression  for  Pk(£,rO  ultimately  reduced  to  a  double  integral  { for  details 
see  C  3 1 ) . 
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SECTION  III 


MULTIPLE  WEAPONS:  INDEPENDENT  PASSES 

1.  Introduction 

It  is  conceivable  that  when  delivering  weapons  to  a  given  target,  an 
aircraft  makes  several  independent  passes  over  the  target  and  releases  one  or 
more  weapons  at  each  pass.  Multiple  weapons  are  assumed  to  be  released 
according  to  a  stick  delivery  pattern.  In  general,  the  number  of  passes  is  a 
random  variable  and  the  number  of  weapons  delivered  at  each  pass  is  a  random 
variable.  The  determination  of  the  probability  of  kill  of  the  target  will 
then  depend  on  the  probability  of  kill  of  the  target  at  each  delivery  as  well 
as  the  statistical  characteristics  of  the  number  of  weapons  released  at  each 
pass  and  the  number  of  passes. 

To  solve  this  problem,  it  Is  necessary  to  know  the  following: 
i.  The  relationship  between  the  probability  of  kill  of  the  target  and  the 
number  of  weapons  released  at  a  given  pass.  This  type  of  functional 
relationship  would  be  the  result  of  an  analysis  of  the  stick  delivery 
of  weapons. 

ii.  The  probability  distribution  function  of  the  number  of  weapons 
released. 

111.  The  probability  distribution  function  of  the  mrober  of  passes. 

In  general,  let 

M  ■  random  variable  denoting  the  number  of  passes; 

N^  «  random  variable  denoting  the  number  of  weapons  released  at  the  ith 
pass  { 1*1 ,2,... ,M) ; 

Pk( Ni )  *  probability  of  kill  for  Nj  stick  delivered  weapons  per  pass; 
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The  net  target  probability  of  kill  is 

pjM  =  1  -  i"1  U  -  W1  (22) 

(Ml 

Even  when  M  and  are  deterministic,  explicit  expressions  for  Pkf'  cannot.  In 
general,  be  obtained  unless  the  function  P^(*)  Is  known.  However,  for  single 
weapon  release  (N,-=l  for  all  i)  and  M  fixed,  a  closed  form  expression  for  Pkf 
can  be  obtained.  This  particular  case  is  investigated  in  detail  in  the 
paragraphs  that  follows'. 

2.  The  Model 

The  development  of  the  model  for  single  weapon  delivery  per  pass  Is 
rather  straightforward  and  consists  essentially  in  determining  the  probability 
of  kill  of  the  target  per  pass  and  then  applying  formula  (22). 

It  wil  be  assumed  that  the  delivery  errors  in  each  of  the  x( range) 
direction  and  y  (deflection)  direction  are  independently  and  normally 
distributed  with  respective  Standard  deviations  o  and  a  .  Let 

a  Jr 

(u,v)  =  cord i nates  of  the  point  target; 

(x,y)  =  cordinates  of  the  point  at  which  weapon  impacts. 

It  is  assumed  that  the  weapon  is  aimed  at  the  target.  Then,  the  probability 
density  function  of  the  delivery  error  In  the  x  direction  is 

f  (x-u) i— «p  (V) 

and  the  probability  density  function  of  the  delivery  error  in  the  y  direction 
is 


f-(y-v)  « '  -  1 . *  exp( -  ■]  (24) 

2  /2Toy  2oJ 
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If  we  assume  that  the  damage  function  is  represented  by  the  three 
parameter  Carleton  damage  function,  then  the  probability  of  kill  at  (u,v) 
given  that  the  weapon  impacts  at  (x,y)  is 

D(u-x,  v-y)  *  Dq  exp{-  DQ[( . ]}  (25) 

Rx  Ry 

To  determine  the  probability  of  kill  of  the  target  for  a  single  weapon  in 
the  presence  of  delivery  errors,  the  laws  of  conditional  probabilities  are 
used  and  one  has 

P^f  *  Probability  of  kill  at  (u,v) 

mm 

3  /  /  [Probability  of  kill  at  (u,v)|weapon  impacts  at  (x,y)] 

[Probability  weapon  impacts  between  (x,y)  and  (x+dx.  y+dy)] 

Using  (23),  (24),  and  (25),  one  obtains 

m  m 

pkf  *  /  J  D(u-x,  v-y)  fj(x-u)  f 2(y-v)  dx  dy 
'  -/*  /*  00  «pi- 

R  R 

*  ,  y 

•  ■£-  dx  dy  (26) 

This  double  integral  can  be  easily  evaluated,  and  it  is  found  that 
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Since  all  passes  are  assumed  to  be  independent,  then  from  (22)  the  net  target 
probability  of  kill  for  M  passes  is 


-  1  -  -  pkf)M 


Using  (27)  in  (28)  results  in 


pW  -  1  -  [1  -  0f 


'Rx  +  2D0  'l 


/Ry  +  2Df)  ‘y 


3.  Estimation  of  E[pj^]  and  Var[P^^j 

The  Taylor's  series  estimation  procedure  will  be  used  to  obtain 

approximate  expressions  for  E[P,^]  and  Varfpj^]  In  terms  of  the  means  and 

variances  of  the  five  input  parameters  Dn,  R  ,  R  ,  a,  and  a  . 

u  x  y  x  y 

Let  nQ ,  Rx,  Ry,  axv  and  7  refer,  respectively, to  the  mean  of  Dq,  Rx,  Ry, 
ax  and  Oy.  Expanding  P^  about  the  point  (O^.  R^,  Ry,  7%,  o*y)  In  Taylor's 
series  and  retaining  only  first  order  terms  one  obtains 


p[J^  (0-,R,R,a,o)*  p[P  (IT- ,  If  ,  If  ,  cT  ,  cT  ) 
kf  .  O’  x’  y*  x*  y  kf  '  O’  x  ’  y*  x*  y 


3p(M) 

3Pkf 


,p(M) 

3Pkf 


+  (D0  -  °0}  30^""  +  (Rx  *  Rx>  nr^~ *  {Ry  *  V 


3P(M)  3P{M) 

/  -»  3Pkf  •  ,  -  ,  3  kf 

,x  x'  3a  v  y  y’  3o 


All  the  partial  derivatives  of  P^  with  respect  to  the  five  variables  Oq,  Rx, 

R„,  a  and  a  are  to  be  evaluated  at  the  mean  values  of  the  variables, 
j  x  y 

/ u\  |u\ 

Approximate  expressions  for  EKf]  and  Var(  P^  J  are  now  provided. 


r 


a 


Estimation  of  E[pj^] 


Taking  expectation  on  both  sided  of  (30)  yields  as  a  first  approximation 


EtPkf}  (0Q,  Rx,  Ry,  V  ay)]  =  *y.  °x»  V 


from  (23) 


E[p[5}]  =  [1  -  Dn 


R*  +  20q  a‘  /  R‘  +  2Dq  op 


For  simplicity  the  bar  (-)  notation  will  be  omitted  in  the  sequel.  Let 


2  2  2 
Q1  “  Rx  +  2n0ax 


2  2  2 
Q2  -  Ry  +  200  °y 


Then,  from  (31) 


1  ,,-In  M 


EC-kf  ]  ■  1  -  »  -  °0  *x  "y  »r  °2  > 


b.  Estimation  of  Var[p£"  ] 

First,  expression  (30)  is  written  as 

Pkf}  {D0‘  Rx*  Ry*  V  ay^  *  Pkf'  V  ®y* 


"W  -  "W  - 

{D0  *  V  3o“  *  (RX  -  Rx}  3r7"  *  {Ry  *  V  iT5~ 


(35) 


,P«? 

*  <»,  *  »,>  -- 


3o 


3P 


a  ) 

y; 


(M) 

kf 


3a 


Squaring  and  taking  expectations  on  both  sides  of  (35)  yields 


aP(M)  ,P(H) 

Varll-W)  =  Vart  D0]  (5^-) 2  *  Var[Rx]  (-^-)2 

,p(M)  lp(") 

*varll.  ],(^f-)2*var[,x] 

y  x 

3P{M) 

9  r  It  f  7 

+  Var[a  ]  (-r- — •)  +  covariance  terms  (36) 

y  y 


Again,  it  should  be  noted  that  all  partial  derivatives  are  evaluated  at  the 
point  (Dq,  Rx,  Ry,  ax,  ay). 

Now,  in  general,  from  (28) 


3P{M) 

3Pkf 

WT 


Ki 


iM-l  3Pkf 


kf 


a  (* ) 


(37) 


Pkf  ^er®  1s  Probability  of  kill  of  a  single  weapon 


Expression  for  3P^/3(«)  can  be  obtained  from  the  results  in  Section  V  in 
which  biased  delivery  error  for  single  weapon  is  considered.  For  the  present 
Section,  results  for  3 (• )  are  obtained  by  setting  In  all  the  partial 
derivatives  of  Section  V  8x*0*By.  This  Immediately  yields  from  (115)  and 
(116)  E1*1»E2. 


From  (117) 


3P 

To: 


kf 


\  j  (1 


On  Or2  *  On  «2  On2] 
0  x  yl  0  y  w2  1 


M 


(38) 


c.  Numerical  Example 

An  aircraft  makes  M=2  passes  over  a  given  target  and  at  each  pass  drops  a 
general  purpose  bomb.  The  weapon  Is  aimed  at  the  target.  The  parameters  of 
the  Carleton  damage  function  are 

00  =1,  Rx  =  120  ft  and  Ry  -  200  ft 

These  parameters  are  subject  to  estimation  errors  and  their  variances  are 
given  by 

Var[D()]  *  .002,  Var[Rx]  =  50  ft^.  and  Var[Ry]  3  70  ft^. 

The  aiming  errors  in  the  x  (range)  and  y  (deflection)  directions  are 
independently  and  normally  distributed  with 

ax  =  30  ft  ;  <Jy  =  20  ft. 

Var[a  ]  ■  4  ft VarCo  ]  «  4  ft2 

*  y 

we  have 

.  Q2  =  (120)? '♦  (2)(1)(30)?  «  16,200 
Qj  '■  127.279,22 

(200)2>  (?)<  1)(20) 2  *  40,800 
02  *  201.990,10 

V  ■  °o  ",  ",  It*'  9j* 
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w-  V 


...  .  jlH12g)(200) 

(1 2?. 2/9, 22) (201.990,10) 


.93352 


From  (34) 

E(P<f]  -  1  -  (1  -  Pkf)2 

*  1  -  (1  -  .93352) ^  =  .9956 


From  (38) 

3Pkf  _  (120 

3^“  ~  (127.279,22 


)(200)  ...  (1  _  (1)(30)2  _  (1)(20)21 

)( 20 1.990, 10)  }l  16,000  40,800  J 


*  .933,52  [1  -  .055,555,5  -  .009,803,8] 
-  .872,505,7 


From  (39) 

I!*! , . . -D.i(2oo) _ r  fl .  jimi] 

( 1 27 . 27  9 , 222 )  ( 20 1 .990 ,  10 )  1 1  16,200  * 

»  .007,779,3  (1  -  .888,888,8) 

«  .000,864,3 


From  (40) 


3Pkf  ( 1)  f  120) _ r.  (200)2i 

(127. 279, 22) (201. 990, 10)  **  *  40,800* 

«  .004,667,6  (1  -  .980,392,1} 

«  .000,091,5 
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( 127.279,22) (201. 990, 10) 3 
=  -  .000,915,2 


Assuming  all  the  covariance  terms  to  be  zero,  then  using  (43)  we  obtain 
Var[P*2)]  =  2  (1  -  .93352)  { (.002)  (.872, 505, 7)2  +  (50)(  .000.864.3)2 


SECTION  IV 


LASER  GUIDED  BOMBS 


1.  Introduction 

For  laser-guided  bombs,  the  distances  between  the  mean  points  of  impact 
( MP I )  and  the  actual  impact  points  cannot  be  accurately  represented  by  the 
circular  or  Rayleigh  distribution  function.  It  becomes  necessary  to  devise 
some  other  distribution  function  to  represent  the  delivery  errors  distribution 
in  the  range  and  deflection  directions.  Let: 

ox  =  standard  deviation  of  delivery  errors  in  the  x  or  range  direction; 

<Xy  *  standard  deviation  of  delivery  errors  in  the  y  or  deflection 
direction; 

FR(r)  *  the  distribution  function  of  the  delivery  errors  R  in  the  radial 
direction; 

a  =  a  positive  constant  such  that  0  <  a  <  1. 

Then,  it  has  been  found  experimentally  [1]  that  FR(r)  is  more  accurately 
represented  by  a  function  of  the  form 

2  2 

F„(r)  *  a  [1  -  exp(-  -£=-)]  +  (l-«)  (1‘  -  exp(-  ~x)]  ,  0  <  r  <  •  (44) 

K  .  ’  2*y  ' 

FR(r)  is  a  proper  distribution  function  since  F(0)  *  0  and  F(«)  •  1. 

The  corresponding  density  function  is 

2  2 

fR(r)  *  exp(-  — £*•}  ♦  eXp(-  -Ay-) ,  0  <  r  <  •  (45) 

°x  2<j  «y  2oy 

In  order  to  find  the  joint  density  function  of  the  delivery  errors  X  and 
Y  in  the  range  and  deflection  directions,  say  fx.Y^***)*  **  *5  necessary  to 
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make  certain  assumptions  about  the  random  variables  R  and  9  which  represent 
respectively  the  delivery  errors  in  the  radial  direction  R  and  the  delivery 
errors  in  the  argument  9.  The  following  shall  be  assumed: 

a.  The  delivery  error  is  equally  likely  to  be  at  any  point  in  the  interval 

0  <  0  <  so  that  the  marginal  density  function  of  9  is 

f9  3  27*  0  <  9  <  2»  (46) 

b.  The  random  variables  R  and  9  are  independently  distributed. 

The  joint  probability  function  of  R  and  9  is: 


fR(r)  f0(0) 


=  *  exp(‘  L2 1  *  IxIr  "p('  J2>1  • 


0  <  r  <  •,  0  <  9  <  2* 


(47) 


The  probabil ity  that  the  weapon  will  impact  In  the  interval  (r,  r+dr)  and 
(9 ,  ,0+d0)  is: 


■  Vr)  dr  de  *  27  [-f-  exp(-  -~) 


2o 


+  llyii  exp( — £*-)]  h  dr  d9 ,  0  <  r  <  -  ,  0  <  0  <  2* 
o'  2o„ 


(48) 


Changing  to  Cartesian  coordinates,  one  obtains  the  probability  that  the  weapon 
will  Impact  in  the  interval  (x,  x+dx),  (y,  y+dy)  or 

'  2_2, 


fx>y(*,y)  dx  dy  *  {-§-  exp{-  ^ 


2"  V 

x 


.2^2 


2<? 


♦  exp[  -  (*■  ^  ] }  dxdy,-»<x<»,-»<y<» 


2a 


(49) 
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Thus,  the  joint  probability  density  function  of  X  and  Y  is 


fx  w(x.y)  -  «  —  exp[ -(-—■)]  •  exp[-(-*jp)j} 

Z2.  ax  2a*  /Zt  ax 


+  (l-«)  {  _1--  exp[^(-~-)]  •  - expK-^y-)]} 

SE  ay  20y  ^2*  ay  2a  y 


0  <  a  <  1,  —  <  x  <•,■—<  y  <  - 


(50) 


Clearly  X  and  Y  are  not  in  general  independent.  X  and  Y  are  independent  only 
in  the  special  cases  when  <x«0  and  o*l.  However,  fy^x.y)  is  a  proper  density 
function.  It  may  be  verified  that  fx.Y***^  *s  a  un "* niod a  1  function,  achieving 
a  maximum  at  the  point  (0,0).  It  may  be  noted  that  ^x.Y^**^  is  9iven  as  a 
convex  combination  of  density  functions  of  two  pairs  of  independently 
distributed  random  variables. 

Let  the  MPI  of  the  weapon  be  the  coordinates  (u,v)  of  the  target,  in 
other  words  the  weapon  Is  aimed  at  the  target.  Let  (dx,  dy)  be  the 
infinitesimal  rectangle  close  to  the  point  (x,y)  at  vtolch  the  weapon 
Impacts.  Define  the  random  variables  X  and  Y  which  measure  respectively  the 
distances  between  the  target  pol.it  and  the  weapon  impact  point  along  the 
abscissa  and  the  ordinate.  Then,  from  (50)  the  probability  that  the  weapon 
will  Impact  In  the  rectangle  d<  dy  Is  given  by: 

fx  Y(x-u,y-v)dxdy  *  a  exp(-  (x‘u? '  * 

’  /2t  a‘  2dx 

exp[-  (x-u)*  +  Util-]}  dx  dy  (51) 

2®. 


♦  (I-*)  {-S' 
2* 


2.  The  Model 

To  determine  the  probability  of  kill  of  the  target  for  a  single  weapon  in 
the  presence  of  delivery  error,  th  laws  of  conditional  probabilities  are  used 
and  we  have  the  probability  of  kill  at  (u,v)  as 


3  j  /  [Probability  of  kill  at  (u,v)  |  weapon  impacts  at  (x,y)] 
[Probability  that  the  weapon  impacts  between  (x,y)  and 
( x+dx,  y+dy)1 

»  cm 

3  /  J  D  (u-x,  v-y)  fx>y(x-u,  v-y)  dx  dy 


=  /  /  0Q  exp{-  DqK-^)2  ♦  (-^)2] 

XV 


(®  1 - ?  ®xp[  - 

2n  a- 

X 


:x-u)2  + 


+  (1-a)  { - -  exp[-  ^ x~u^ — } )  dx  dy 

2n  o'  2<jc 

y  y 

Making  the  change  in  variables  w  »  x-u  and  z  »  y-v  yields 


2  2 
w  .  z 


2  2, 
w _ »  z) 


r  n0  expt-  ♦  -^j}  (a  l-^-  exp{-  . -')} 


-  -  -  Rx  «;  2-  °x 

?  2 

*  (l-a)  { - exp[--^ — -  *-^]})  dx  dy 

2*  °y 


This  integral  may  be  expressed  as  a  sum  Involving  products  of  two  single 
integrals: 


«  °n  *  °n  "  J 

pkf  “  ~ expt^“TT  +  7T 

2*°x  -  Rx  2ox 


V \  z2 


♦  T-)]  dw«  /  *xp{- (“V*  +  “V'))dz 

2a:  -  K,  2a: 


y  •  x 
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(1-a)  0 


0  ■ 


2*o 


°n  w  J 

e*p{  -  (““2  +  ~JV‘)1  dw*  / 

Rx  2<Jy 


exp(-  +  -£-j)]  dz 

«y  *} 


.2 


Consider  any  one  of  the  four  Integrals.  For  example,  let 


(54) 


■** '  I  “4  -  C-f  ♦  rf)  «21  d« 

-  Rx  2ox 


n. 


l 


Let  *  (-2.  +  )2  w 

/7.  Rx  2°x 


(5,5) 


(56) 


Then 


',X  -  -pr  ♦  "V*  /*« P(-4* 

xx  /2  Rx  2<*x  0  ‘ 

.->-A  /ST 

.  /7  R  x  2®x 

_1 

■  ^  «,  V  (<  *  °2S* 


(57) 


The  other  three  Integrals  may  be  obtained  In  a  similar  fashion.  One  then 
immediately  obtains  for  P^ 


\t  •  •  °0 


J<  *  *°0  »x  /rJ  *  2tl0  »x 


♦  (1-.)  0 


®  n— 

/R"  ♦ 


20„  a2  /r2‘  ♦  20-  a2 
x  0  y  y  o  y 


3.  Estimation  of  E[Pkfl  and  Var[P^} 


(58) 


Recall  that  P^f  Is  a  function  of  the  six  parameters  a,  0Q,  R^,  Ry,  <?x  and 
ay.  let  a,  00.  Rx,  Ry,  o%  and  refer  respectively  to  the  mean  of  a,  HQ., 
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Ry.  °x  and  a  .■  Expanding  P^  about  the  point  (o,  Dg,  Rx,  Ry,  ox  ,  «y)  one  obtains 
up  to  the  first  order  terns: 

pkf  (“•  °0*  Rx‘  Ry*  °x»  V  *  P*f  *°*  ^0*  *x*  Ry*  °x* 

^Pkf  —  —  ®Pkf 

+  (o^  —  +  (00  -  D0)  -Jo-  ♦  (Rx  -  Rx)  3l^ 

3P.f  _  3P.f  _  3P  - 

*  !Ry  -  v  IT  *  <°x  -»,)-#*  (•,  -  air~  159) 

y  x  y 

Note  that  the  partial  derivatives  of  P^f  with  respect  to  the  six  variables 
a,  Og,  Rx,  Ry,  ax  ,  and  ay  are  to  be  evaluated  at  the  mean  values  of  the 
variables.  We  now  provide  approximate  expressions  for  EtP^f]  and  VarCPj^]. 

a.  Estimation  of  E[Pkfl 

Taking  expectations  on  both  sides  of  (59)  yields  as  a  first  approximation 


E(  Pkf(°*  V  Rx*  Ry*  °x  °y)l  *  pkf(®»  ^0»  *«•  *y»  °x*  °y) 


=  a  D„ 


J  K *  2Do  4  1  "5  *  2Vx 


*  * 


I  ¥y  ♦  2n0  7* 


(60) 


b.  Estimation  of  VarfP^] 

First,  expression  (59)  is  written  as 

pkf(°*  D0»  Rx’  Ry»  ®x»  ®y)  "  Pkf(®»  ^0»  *x*  V  ®x*  ®y) 
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3Pkf  -  3Pkf  -  3Pkf 

<— )  IT  *  (°0  -  V  Tof  *  <»x  '  "x>  3*7 


-P.f  _  3P.f  _  3Pkf 

♦  (Ry  -  Ry)  fR  +  *  °x)  3o  +  (°y  “  °y^  3o  “ 


Squaring  and  taking  expectations  on  both  sides  of  (61)  yields 


j  p  3  p  3  p 

Var(Pkf]  »  varta]  (- -)2  +  VarCDn]  (j^)2  +  Var[Rx]  (^)2 

0  x 


♦  Var[R  ]  (jfl— )  +  Var£ax3  (jo — )  +  VarCoy]  (3^  ’)' 
y  x  y 


+  covariance  terms 


Again,  it  should  be  noted  that  all  the  partial  derivatives  are  evaluated  at 
the  point  (a,  Dg,  R^,  Ry,  <TX#  oy) .  We  now  obtain  the  expressions  for  the 
partial  derivatives 


ST1  “  °0  %  ”y  *  2°0  °x)’J  Oj  *  2°0 


-  \  \  %  (»2  ♦  2  b0  .2)  1  (B2  ♦  2D0  «2)  ! 

-  °0  “x  tR»  *  Z00  ax)'5  c “5  ♦*%»«>'* 


•  \  A  tRx  *  2D0  *  2"o  «S>  *1 
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V, 

K~ 

k.* 


*  fl-“)  *x  Sy  «Sx  *  2I,0  «y)"5  (R2  ♦  ?  »„ 


.3  1 

-  °0  °l  ("x  *  2"o  »J)  Oy**  Vy>'2 

-1  3 

-  °0  °y  ("«  *  200  "y)'J  t»v  *  2°0  »y)"2I 


dp 


kf 

X 


1 


1 


l 


3R  ■  ■  °0  Ry  *®0  °x)’J  !(**♦  200  %V? 

_3 

‘  Rx  (Rx  +  2D0  *  f1**)  n0  Ry  (RJ  +  200  *yf? 

f(Rx  +  2D0  °y)  ?  '  Rx  (Rx  +  2DOffy)”2] 


kf  i 

W  *  “  D0  Rx  (RX  +  2D0  ffx)  ?  f(Ry  +  20o  °x)~? 


(64) 


(65) 


>• 


,  I. 


-  #y  <Ry  *  %  ”x)  Z1  *  <l->  »„  ",  0>x  ♦  ®o  «y)'? 


.2,-2 


o  V 


t(#y  *  2D0  "y>  ?  -  Ry  tRy  *  ®0 
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5!f  "  *  Do  Rx  Ry  I'  2Do  «,  («j  *  ®0  «xVJ  Of  *  ®„  °IP 


y  ”  0  v  'x  ’  ‘■"o  VXJ 

-l  .3 

-  2Dq  o  (R2'+  20q  a2)  ?  (R2  +  2Dq  o2/2) 


(66) 


(67) 


r . 

•  t 


3Pkf  .1  3 

3*y  *  (1-a)  n0  Rx  Ry  t-  200  °y  (Ry  *  2D0  *y)  ?  (R*  *  2n0  °y)*? 


'  ?D0  ffy  (Rx  *  200  *y)  ^  (Ry  +  200  °yf?] 


(68) 
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c.  Numerical  Example 

Assume  that  for  a  laser  guided  weapon ,  the  parameters  of  the  Carleton 
damage  function  are 

00  -  1,  Rx  *  120  ft  and  Ry  «  200  ft 

These  parameters  are  subject  to  estimation  errors  and  their  variances  are 
given  by 

Var[D0]  *  .002,  Var[Rxl  =>  50  ft2  and  Var[Ry]  =  70  ft2 

The  guided  weapon  is  aimed  at  a  point  target.  The  aiming  In  the  x(range)  and 
y  (deflection)  directions  have  a  joint  probability  density  function  given  by 
(51)  where 

a  *  .75;  a x  *  30  ft;  oy  *  20  ft 

Var[al  *  .001;  Var[o  ]  *  4  ft2;  Var[o  ]  =*  4  ft2 

*  y 

To  obtain  the  numerical  values  of  the  expressions  for  ECP^I  and  the 
partial  derivatives,  the  following  quantities  are  introduced 


Qj  *  Rx  ♦  20q  <»x  (69) 

»  Ry  ♦  200  oj  (70) 

03  *  Rx  *  200  °y  <7l> 

0j-Ry  +  2°O’x  {72) 
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Substituting  for  the  values  of  the  parameters  yields 


Q2  -  (120)2  ♦  (2)(1)(30)2  *  16,200 


«  127.279,22 

»  (200)?  +  (2)(1)(20)2  -  40,800 


* 


Q?  =  201.990,10 

Q2  =  ( 120) 2  +  (2) (1) (20) 2  *  15,200 
O3  *  123.288,28 

Q2  *  (200) 2  +  (2) (1) (30) 2  *  41,800 
04  =  204.450,48 


Omitting  for  simplicity  the  symbol  for  averages  one  obtains 


E[Pkf]  *  a  Do  Rx  Ry  °11  Q4X  +  ^1*a)  °o  Rx  Ry  <*3  Q2l 


°0  Rx  Ry  ta  ^  041  +  021] 


( 1 ) ( 1 20 ) ( 200 )  [- 


.  _ j25 _ , 

( 1 2  3 . 288 , 28 ) ( 20 1 . 990 ,10 T 


932,649,1 


Now 


3P 


sr1  ■  "o "« R,  lo;1  o'1  -  q-1  q;1] 

,450,48) 


(74) 


"  (T23.288,'28)(201,990,10)^  *  ”*041,452,8 

3P.  f 

3DT- a  °  Rx  Ry [Q? -  Vx  q;1  -  Vx2  Qil  o;3i 


*  (1-a)  [Q"1  Q*1  _  Oq  g”3  g”1  .  g^  ff2  q-1  q-3j 


Rx  %  *?  OJ1  [  1  -  %  <  9?  -  00  .*  Q-?] 


2  „-2 


.2  n”  2l 


.-1  „-l 


Mr-a)  Rx  Ry  Qjl  Oi1  [1  -  o0  *y  O32  -  00  Oy  Q-?] 

■  ^^7 •  *279^22? ^ 204^450 ^48)  f 1  *  ^[lol  *  Hl'iloo  ^ 

+  ( ZO 1 .990^  10 )  ( 1 2 3 ^ 288 , 28 )  1 1  *  ^joo  '  *  Hofloo  *1 

-  (.691, 714, 4) (.922, 913, 3)  +  ( .240, 934, 7)( .963,8802) 

•  .870,6246 


(75) 


3f> 


kf 


■  *  °0  *y  041  t®il  -  «;3]  *  »->  »0  %  a;1  [o;1  -  r3  q;3] 


*  °0  Ry  °i1  °4l  (1  -  R,  0j2l  *  <l-»)  00  Ry  Oj*  Oj‘ll  -  R3  Oj’l  (76) 


1  o“lr ,  „2  _-2i 
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»>  w>  v>  i' 


*  , 

JV-L’VJV.V:  -v*.  ■\.W\  - 


l.A.  '  ^ 

’M-.T  l-’AW. 


. Ljsnmm-.  T r i .  im2.) 

(127. 279. 22)  {204. 450, 48)  11  16,200J 

(.ZS)(l)(gQ0)  ,  (120)2, 

(201.990, 10){123. 288, 28)  l1  15,200* 

-  (.005, 764,2)(. 111,111,1)  +  (.002,007,7)(. 052,631,5) 

*  .000,746,1 

3P 

3r  =  af)o  Rx  ^  [tf-  Ry  q;3]  ♦  (IhO  d0  rx  -  R2  o;3] 

a  °  °0  Rx  <5*  [1  -  Ry  Ol2]  ♦  (Ihx)  00  Rx  Q*1  Q51  [1  -  Ry  Q^l  (77) 

,  _ _ (_.75)(1)(120)  f  (200) 2i 

( 127 .279.22 )  ( 204.450,48)  «■ 1  41 ,800J 

+  .  ( *25) (1) (120)  r  r .  ( 200) 

■  (201.990,10) ( 123. 2o8, 28)  1 1  '  40,800* 


*  (.0003,458, 5) (.043,062,2)  +  ( .001 ,204,6) ( .019,607,8) 
=  .000,172,5 


3“ 


nf  *  - z*  °o  % R*  \  i«;‘  q'i  *  n;1  °;3i 


2a  Og  <JX  Rx  Ry  Q^1  Q*1  [q‘2  +  qj2] 


.  iilLZi )I(1)2(30)(120)(200)|  r . 1 . +  1  | 

(127.279,22) (204.450.48)  L 16,200  41, 800 J 


.003,554,7 


(78) 
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(79) 


3  P 

-  -  2(l-a)  D2  oy  Rx  Ry  (O"1  Qj3  +  Qj1  Qj  ] 


-  -  2(l-a),  Dq  ox  Rx  Ry  Q‘2l  O31  [O32  +  Og2] 

-  -  .000,870,2 

To  determine  the  variance  of  P^f,  we  assume  a,  Og,  Rx,  Ry,  ox  and  ay  to 
be  Independent  so  that 


3  P  3  p  3  P 

VarCPkf]  «  (Varfa])  (-J£)2  ♦  (Vartrig])  (j^1)2  ♦  (Var[Rxl)  (^)2 

0  ^ 


+  (Var[R  ])  (57^  2  +  (VarOx])  (37^) 2  ♦  (Var[o  ])  (37^ 2 
J  y  x  y 


3PL 


*  (.001) (-.041, 452, 8) 2  +  ( .002) ( .870, 624, 6)2 
+  (50)(.000,746,1)2  +  (70)(.000,172,5)2 

♦  (4)(-.003,554,7)2  +  (4)(-.000,870,2)2 
«  .001,601,1 

aQ  m  .0400  J _ 

Pkf 


*  two-standard  deviation  confidence  Interval  on  Pkf  Is 

•  .9326  t  .nsoo 

He  note.  In  particular,  that  the  contribution  to  the  total  variance  of  Og  is 
94.7*. 


SECTION  V 
GUIDED  WEAPONS 


1.  introduction 

In  the  most  recent  JTCG/ME  manual  on  the  "Derivation  of  JMEM/AS  Open-End 

Methods"  (1983)[21,  the  methodology  for  computing  range  and  deflection 

delivery  errors  for  guided  weapons  Is  outlined.  It  Is  assumed  that  In  the 

range  direction,  the  aimpoini  is  offset  by  a  negative  value,  that  is  the 

center  of  the  distribution  of  impact  points  Is  short  of  the  aimpoint.  No  such 

offset  occurs  in  the  deflection  direction.  The  amotnt  of  offset  is  a 

correction  factor  which  can  be  Introduced  as  a  bias  In  the  distribution  of 

delivery  errors  in  the  range  direction.  This  bias  Is  assumed  to  be  a  known 

input  to  the  problem.  The  distribution  of  delivery  errors  in  each  of  the 

x(range)  and  y  (deflection)  directions  are  assumed  to  be  normally  and 

independently  distributed  with  the  same  standard  deviation  (a  .  *  <x  *  a). 

*  y 

The  problem  treated  in  this  section  Is  quite  general  and  considers  the 

presence  of  biases  in  both  the  range  and  deflection  directions  with  »x  * 

8y  setting  the  bias  in  the  deflection  direction  equals  to  zero  and 

at  *  a  =  a,  the  distribution  of  peli/ery  errors  for  guided  weapons  is 
x  y 

obtained.  By  setting  the  biases  in  both  direction  equal  to  zero,  the 
distribution  of  delivery  errors  for  general  purpose  bombs  is  obtained. 


2.  Weapon  Delivery  with  Offset  or  Bia 

In  studying  the  delivery  of  a  guijled  weapon  (idealized  as  a  point)  upon  a 


point  target,  several  factors  have  to  l>e  accounted,  for  in  developing  an 
expression  for  the  probability  of  kill  of  the  target.  Among  such  factors,  the 
coordinates  of  the  various  interplaying  variables  should  be  considered. 
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particularly  if  the  weapon  Is  not  aimed  directly  at  the  target  or  the  presence 
of  bias  In  the  delivery  of  weapon  cannot  be  avoided. 

Let  then; 

(u,v)  »  location  of  the  point  target  on  the  ground  plane; 

(C,n)  =  coordinates  of  the  desired  mean  point  of  impact  (HPI)  of  the  weapon, 
or  the  coordinates  of  the  point  at  which  the  weapon  is  aimed  or 
aimpoint; 

(x,y)  =  coordinates  of  the  actual  weapon  impact  point; 

(i,m)  «  coordinates  of  the  center  of  impact  points  (CIP). 

A  little  elaboration  is  needed  on  the  last  set  of  coordinates  defined. 

If  there  is  no  bias  in  the  delivery  of  the  weapon,  then,  the  CIP  coincides 
with  the  MPI.  If  on  the  other  hand  bias  exists,  then  the  CIP  will  be 
different  from  the  MPI.  Thus,  in  general,  I  *  5  and  m  *  n. 

Assume  now  that  the  x  direction  is  the  direction  of  the  range  and  the  y 
direction  is  the  direction  of  deflection.  Suppose  now  that  there  is  a 
delivery  error  when  releasing  the  weapon.  This  error  is  measured  in  each  of 
the  x  and  y  directions.  The  delivery  error  in  the  x  direction  is  the  distance 
along  the  abscissa  between  the  weapon  impact  point  and  the  center  of  impact 
points  or  (x-i).  Similarly,  the  delivery  error  in  the  y  direction  is  the 
distance  along  the  ordinate  between  the  weapon  impact  point  and  the  center  of 
impact  points  or  (y-m).  Assume  (x-i)  and  (y-m)  to  be  independently  and 
normally  distributed  with  respective  standard  deviations  <»x  and  oy.  Thus, 
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(1M,)  ’  /jr 


1 


exp[ 


h=sl! 

2o? 


•] 


(81) 


y  -  y 

Assume  that  the  bias  in  the  x  direction  and  the  bias  in  the  y  direction  are 
known.  Thus,  let 


ry 


bx  *  5-t  =  bias  in  the  x  direction  or  range  offset, 
by  =  n-m  =  bias  in  the  y  direction  cr  deflection  offset. 


Thus, 


l  =  £-b 

■  X 

m  =  n-b 


(82) 


and  it  becomes  possible  to  express  the  probability  density  functions  of  the 
delivery  errors  in  terms  of  the  MPI  and  the  offset.  Thus 

>2 

( »_r  aK  ^  t*  mi  1 1  ■  avnT  —  .  * 

1 


.  (x-e+b  ) 

f,  (x-e+b  )  »  — : -  exp[- 

x  /&T  a 


2a  ‘ 


■] 


, ,  (y-n+b  ) 

<y_0+t>y)  - 


(83) 

(84) 


3.  The  Model 

If  we  assume  that  the  damage  function  is  represented  by  the  Carleton 
damage  function,  then  the  probability  of  kill  at  (u,v)  given  that  the  weapon 
impacts  at  (x,y)  is 


D(u-x,  v-y)  »  0Q  exp{-  Dq  [(^)2  ♦  (^)2]} 

x  y 


(85) 


To  determine  the  probability  of  kill  of  the  target  In  the  presence  of  delivery 
error,  the  laws  of  conditional  probabll itles  are  used  and  we  have: 
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1 


*  Probability  of  kill  at  (u,v) 


m  m 

-  ^  J. 

«•  «• 

-  /_  J, 

a*  « 

■  J  / 

«a*  -«• 

.  1 
'*°x 


[Probability  of  kill  at  (u,v)|weapon  Impacts  at  (x,y)]. 
[Probability  weapon  Impacts  between  (x,y)  and  (x+dx,  y+dy)]. 
D(u-x,  v-y)  fj  (x-C+bx)  f2<y-n  +  by)  dx  dy 


In  general,  thus,  P^f  Is  a  function  of 

a,  the  location  of  the  target  (u,v); 

bs,  the  coordinates  of  the  MPI  U,n); 

c.  the  biases  bx  and  by  in  the  x  and  y  directions. 

To  evaluate  the  double  Integral  In  the  expression  for  Pkf,  write  P^f  as 
the  product  of  two  single  integrals. 


pkf  *  %  * 


1 


1 


/5T  o%  oy 


•  /  exp{-  2 - + - 7 - ]}  dx 

—  to' 


*  /  exp{-[ 


00(v-y)2  (y-n  ♦  by)? 


2 o 


■]}  <*y 


(87) 


let 

w  *  x-C  +  t>x 
z  «  y-n  +  by 

Then 


jo  r  r«  +  bx)-w)‘ 


pkf *  2^~r  I  expM°0(* 


x  y  — 


■)  *  -=yi  *• 

*x 


.  . .  .(v-n  +  b  )-z  ?  J 

/  exp(-[00( - - f  + 

-*•  y  2o 


■]}  dy 


R  R  ,  u-c  +  b  2  v-n  +  b  „ 

x-£-  «xp(-tl((-^)  ♦C-^)2] 


*  Pkf(u-S  ♦  bx,  v-n  ♦  b  ) 


where 


r2 

2  Kx  A  2 
ql  "  2Dq  *  °x 


?  y  2 
q2  = 


4.  Special  Cases 

Several  special  cases  are  now  considered 
a.  Case  of  no  bias 

Here  bx»0  and  b*0  and  one  obtains 


v. 


(88) 

(89) 


(90) 


(91) 


(92) 


(93) 
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b.  Case  of  no  bias  in  the  y  direction  (guided  weapons) 

This  Is  the  case  corresponding  to  guided  weapons  and  the  expression  for 
Pkf  becomes 


RvRu  i  u-c  +  bw  9  ~ 

pk,  C-MV  v-n)  -  Mp(.  I  [(-^-JS.)2  ♦  (^)2]| 


n 


c.  Case  of  no  bias  in  the  y  direction  and  weapon  is  aimed  at  target. 
Here,  the  desired  MPI  coincides  with  the  point  target  (u=£,  v-n): 

»kf  <»«.  by)  -  ^  «Pt-  1  [(J.)2  ♦  (^)2]J 

5.  Estimation  of  E[Pkf]  and  VarfP^^] 


Let 


\9  u<  *  \ 

*  v-n  +  by 


The  expression  for  P^  becomes,  in  general. 


V  ‘3x- 


vrfiere 


2  x  ? 
qi  ■  roj +  9x 


4  •  rfc 


♦  o 
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(94) 


(95) 


(96) 

(97) 


(98) 


(99) 


(100) 


AV- 

V 

V 

S.  -l 
n 

Iv; 


a u 

4** 


£ 


£ 


fe: 


-%  V 


& 


Note  here  that  Pkf  Is  a  function  of  seven  parameters  namely  Dg,  Ry,  Ry,  o^, 

0^,  Bx  and  By.  We  shall  assume  that  8X  and  By  are  precisely  known  and  are  not 
subject  to  estimation  error,  but  that  the  five  other  parameters  Dg,  Rx,  Ry, 

°x  and  °y  may  be  subJect  to  estimation  error.  It  shall  be  assumed  that  both 
the  means  and  the  variances  of  the  parameters  are  known.  First  we  obtain  an 
explicit  expression  for  P^.  Substituting  for  q^  and  In  P^  (Bx,  By)  yields 


f  I 


\ f  (*,.  V  - 


x  2 
2  00  +  °x 


2  °o  *  °y 


-  exp  {-  \  (■ 


4-+a2 
°0  x 


2  On  +  % 


*1  * 7  vl 


exp  l-\  (-j 


2  D0»2 


R;  * 2  v; 


Ry  2  D0  °y 


exPl-  l  ( 4  ^  p] 


(101) 


The  Taylor’s  series  estimation  procedure  will  be  used  to  obtain  approximate 

expressions  for  E£Pkf]  and  Var[Pk^3  In  terms  of  the  means  and  variances  of  the 

five  Input  parameters  Dg,  Ry,  Ry,  a  and  o  . 

y  x  y 

Let  n0,  Rx,  Ry,  7X  and  7y  refer  respectively  to  the  mean  of  Og,  Rx,  Ry, 

%  and  <y  Expanding  Pkf  about  the  point  {Og,  Rx,  Ry,  «  ■)  one  obtains  up 
to  the  first  order  terms: 


Pkf  (00*  Rx*  Ry*  V  V  "  Fkf  <V  *x»  *y»  V 
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(102) 


+  (o0  -  o0)  3D 


+  <Rx“Rx>«7+<Ry-  Vaftf 


-  3Pkf  -  aPkf 

*  y 


Note  that  the  partial  derivative  of  P^  with  respect  to  the  five  variables  Dg, 
Rx.  Ry,  ox  and  are  to  be  evaluated  at  the  mean  values  of  the  variables. 
Approximate  expressions  for  ECP^f]  and  VarCP^f]  are  now  provided. 

a.  Estimation  Of  ECP^] 


Taking  expectations  on  both  sides  of  (102)  yields  as  a  first  approximation 


E[Pkf(®0*  Rx»  Py*  ax*  ®y)l  *  ^kf  (Vo.  *x-  ^y»  °x*  °y) 

1  B  A  By  A 

*  rr*r  «p{- \  [(^-)2  +  (^-)2]}  (103) 

2ql«2  qi  q2 

(104) 

(105) 

First  expression  (102)  Is  written  as 

pkf  ^n0*  Rx*  Ry*  V  *  pkf  **•  ®x*  °y^ 

3Pkf  _  3  Pfcf  3P., 

“  (00  *  V  30^"  *  {Rx  '  Rx>  3R^  ♦  <Ry  *  fy 
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where 


R2 

-2  Kx  —2 

\-Tir0*°\ 

-2  S_  -2 


b.  Estimation  of  Var[Pkf] 


(106) 


3Ptf  _  aPkf 

'♦  K  -ax^  3 r~+  (oy  "V  XT' 


Squaring  and  taking  expectations  on  both  sides  of  (106)  yields 


3P^f  *  ? 

kJj\Z  ,  <>.  .rn  i  f _ 


*"tPkf]  -  Var[D„]  (jjHj'  *  V.r[Rx]  (j^-) 

U  A 


r  b 


3P^r  9  ? 

+,  Var[R  ]  (gj— )  +  Var[oxl  (gj— ) 


+  Var[o  ]  — ■) *  +  covariance  terms 

y  y 


(107) 


Again,  it  should  be  noted  that  all  partial  derivatives  are  evaluated  at  the 

point  (1L,  IT  ,  R"  ,  <T  ,  of  ) .  We  now  otbalrt  the  expressions  for  the  partial 
u  x  y  x  y 

derivatives 


a.  aV/,00 


ilif  •  lRx  l’,2  *  2  °0  «2’'5  «rf-  °0  Bx  (Rx  *  2  "o  “x)'1! 


l«y  *  2  °0  °y)"2  exp{.  O0  B2  ( R2  ♦  2  t)Q  o*)A\\ 


+  {-  °0  Rx  .*  (Rx2  *  *  D0  2  exp  l-  °n  (R«  +  2  Dn  *5 ‘*1 


0  x  ^x  c  w0  V 


Ry  (Ry  *  2  n0  aP  1  expt*  D0  BbRy  +  2  °0  0y}"ll} 


{f>3  Rx  (Rx  +  2  °0  °x)  ^  I*  Bx  (Rx2  +  2  °0  °x)'1  +  2  °Q  8x  tfx  (Rx  +  2  tyf)”*] 


•  exp(-  0Q  B 2  (R2  ♦  2  00  a2)'1]  •  Ry(R2  ♦  2  D0  2 
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*  U«y  ♦  2  0Q  op  ‘  exp[-  00  B2  (R2  +  2  0Q  op'1] 


Ry  <Ry  +  2  °0  °?  *  extf-  D0  By  (Ry*  +  2  °0  ^ 


2  fn2 


a  n2  f n2 


2  frZ 


+  2  Do B; R;  (r;  +  2  y;)  o$  *  2  d0  ap 


r.  v, 

i r  V 


f  fe.. 


*  e*P(-  D0  B2  (R2  +  2  D0  a2)-11 


(110) 


d.  3Pkf/3°, 


3o~  *  °0  Rx  Ry  (Ry  +  2  Pq  ay)  *  «P(-  D0  By  (Ry  +  2  °0  0y)'^ 


y  v  y  w  “0  yJ 


•  (  -?  D0  ax  (R2  ♦  2  0„  ,2)  ?  e«p(.  0„  B2  (R2  ♦  2  Dg  a2)'1] 


1 

+  4  D0  Bx  ax  (Rx  +  2  °0  °x^"2  (Rx  +  2  °0  ax^’2 


expf-  °n  BX  (Rx2  +  2  °0  *x)'1]} 


(111) 


3 P  73o 

e.  '  icr  j 


5JT  *  Do ",  ",  ("x  *  2  ”0  •&  ?  “d-  °0  BX  (",  +  2  °0 


.3 

•  I-2  D0  ",  *  »0  ’p  ’  n0  *1  ("y  *  2  °0  "J)'1] 


+  4  °0  RJ  °y  (Ry  +  2  °0  "y)‘?  2  T30  ,2)"2 


exPl-  %  By  (Ry  *  2  °oap“H 


(112) 
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/  . 


c.  Numerical  Example 

Assume  that  for  a  guided  weapon  the  parameters  of  the  Carleton  damage 
function  are 

Og  -  1,  Rx  *  120  ft  and  Ry  *  200  ft 

These  parameters  are  subject  to  estimation  errors  and  their  variances  are  given 
by 


Var[Dg]  *  .002,  Var[Rx]  =  50  ft2  and  Var[Ry]  *  70  ft2 

The  guided  weapon  is  aimed  at  a  point  target.  Thus, 

u  »  $  and  v  *  n 

The  aiming  errors  in  the  x  (range)  and  y  (deflection)  directions  are  normally 
and  independently  distributed  with  the  same  standard  deviation 

a  "  a  *  a  *  20  ft 
x  y 

and  Var[o]  «  4,  ft2 

The  bias  in  the  range  direction  is  bx  «  20  ft.  There  is  no  bias  in  the 
deflection  direction:  by  «  0  ft.  It,  thus,  follows  that 

Bx  *  u-C+bx  «  20  ft 

8y  »  v-n+by  «  0  ft. 
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Assume  that  all  the  covariance  terms  are  zero.  For  simplicity,  the  bar  (— ) 
notation  for  expectation  will  be  omitted. 

To  obtain  the  numerical  values  for  and  the  partial  derivatives,  the 

following  quantities  are  introduced: 

Q?  -  Rx  +  2  °0  ax  (H3) 

Q2  =  Ry  +  2  D0  ay  UM) 

EiaexP{-00  BX  [Rx  ^  2  °a  ,  . 

■  exp[ -  Dg  B2  Qj2]  (115) 

E2  -  exp{-  D0  82  [R2  +  2  Og  a2]'1} 

»  exp[  -  00  82  Q’Z]  (116) 

Substituting  for  the  values  of  the  parameters  yields 

Q2  -  ( 120 ) 2  ♦  (2)(1)(20)2 

«  15,200 

Thus,  Qj  *  123,288,28 

*  ( 200) 2  ♦  (2)(1)(20)2 

«  40,800  .  .. 
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Thus,  Q2  «  201.990,10 


Ei '  e,pt  -  1Hfi§o  1  *  -”4-027-4 

E,  .  expf ,  ,U)(0)21  .  i 
2  e  Pl  40,800  J  1 


Using  (113),  (114),  (115)  and  (16)  In  (103),  the  following  expression  for 
ECPfcf]  Is  obtained 


'  “o  "x  »l‘  E1  Ry  <  h 

.  (1) (120) (.974 ,027, 4) (200) (1) 
( 123. 288, 28) (201. 990, 10) 

*  .938,708 


The  partial  derivatives  are  now  computed.  They  are  first  rewritten  using 
(113),  (114),  (115)  and  (116). 

From  (108) 

is*  *  K  ei  «y  »;*■  y  - 1- »» "x  °l  oj3  ei  %  '’j1  y 

+  ((-  Dq  Rx  B2  Q'3  +  2  Dg  B2  Ry  a2  Q“5]  Ej  R  Q*1  E2 


*  »-  «o  «*  ci-  -J  ea>  , 

*  l°0  "x  ®il  EV  «y  Qil[ -By  Qj2  ♦  2  D0  B,  »2  Q|4]  E2I 

-  *x  R  QJ1  QJ1  Er  E2  { 1-  Dg  o\  Qj2  -  0Q  B2  Qj2  ♦  2  O2  82  «2  Qj4 


-  Dq  op  Op  *  D0  Bp  0 p  [-1  +  2  0 Q  op  Q p]\ 


3  (120)(200)(.974,027,4)(ll  , 
JUp  ( 123.288,28)  (201.990,10)  1 

*  (2)(1)2(20)2(20)2  ( 1 ) ( 20 ) 2i 

»  ■  40,800  3 


(1)(20)2  (1)(20)2 
15,200  '  15,200 


(15,200)' 


=  .938,708  [  1  -  .026,315,7  -  .026,315,7  +  .001,385  -  .0f;i,803,9] 


=  .881,399,5 


From  (109) 


■  no  %  021  s  i0!1  ei  -  *«  o;J  h* 2  %  <  k  «r  Ei! 


2  ,>-3 


,2  .2  rt-5 


3R.  0  y 


nO  Ry  <>2  E1  =2  -  °r  <  J  °0  Bx  Opl 


-2  .2 


,2  „-2-> 


'  1)  (200)  (.974.027 .4)  (1)  r,  (120P  (2)(1)(20 

!  1 23 . 288,28) ( 201.990 , 10 )  1  1  *  15,200  **  15,200 


.007,922,5  [ 1-  .947,368,4  (1  -  .052,631*5)] 


.000,801,7 


From  (110) 


,2  02 


XT  ■  “o  Rx  <  E1  1°?*  E*  -  *1  .<  E2  *  2  °0  "y  °j£  M 


°0  Rx  Qll  °2l  E1  E2  u  -  Q2'  Ry  U  -  2  %  K  <£*)] 


-2  „2 


,2  .-2* 


3 


.004,693,5  (1  -  .980,392,1) 


.000,092 


From  (111) 


dp 

5sf  ■  »o  Rx  R,  °1  E2  °0  °x  o;3  E1  ♦  4  °0  8x  ",  "I5  V 


■  2  "o  Rx  %.  °i3  °21  °x  E1  ri  (->  *  2  “o  8x  °12) 

.  (2)(1)2(120)(200)(20)(. 974.027 .4)(1)  ,  j  |(2)(1)(20)2i 
(123.288,28) 3  (201.990,10)  15,200 

-  .002,470,2  (-1  ♦  .052,631,5) 


-.002,340,1 


From  (112) 


3p 

STf  -  D0 >x  ",  "I1  E1  <-2  °0  ",  »23  E2  ‘  *  "S' »25  E2> 


2  °0  R,  R,  E!  E2  ( ->  *  2  "o  By 


.  l2Hl)2(I20)(200H20V(.W.fl27.«)(I) 
( 123. 288, 28) (201. 990, 10) 3 

•  -.000,920,3 


(120) 


(121) 
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To  determine  the  variance  of  P^f,  expression  (107)  is  used 

3  P  3  P  3p 

Var[Pkf]  =  VarCOg]  (g^)2  +  VarCRx]  (g^)2  +  Var[Ry]  (jiT^2 

0  x  y 

*  V*r[o  ]  (^Si)2  *  V«rC.  3  (^*i)2  (12«  «• 

x  y 

*  (.002)(.881,399,5)2  +  (50)( . 000,801,7)*  ' 

*  ( 70) ( .000.092)2  +  (4)(-.002,3401)2 
+  (4)(-.000,9203)2 

=  .001,611,7 

Thus,  aD  *  .0402 

Kf 

A  two- standard  deviation  confidence  Interval  on  E[P|(f]  is 

Pkf  *  E^Pkf^  *  2  aPkf 
*  .9387  i  .0804 

a 

It  is  noted,  in  particular,  that  the  contribution  of  VarfDgl  to  the  total 

variance  Is  about  96. 4t.  1 


SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  report,  the  modeling  of  multiple  weapons  delivery  has  been 

attempted.  Details  related  to  the  stick  delivery  of  weapons  are  described  In  a 
A 

separate  report  [3],  In  addition,  the  modeling  of  the  delivery  of  guided 
weapons  has  been  undertaken  under  two  differnt  sets  of  assumptions  for  the 
distribution  of  delivery  error. 

Although  the  underlying  methodology  for  tackling  these  different 
situations  is  the  same,  nevertheless,  the  statistical  characteristics  of  the 
distribution  of  delivery  errors  constitute  the  main  element  which  brings  forth 
the  differences  between  the  various  models. 

The  solutions  to  the  problems  related  to  multiple  weapons  are  far  from 
being  complete.  For  example,  the  approximate  solution  to  the  general  stick 
delivery  problem  as  partially  reported  in  [3]  needs  some  verification  using  the 
exact  approach  presented  In  Section  II  In  this  report.  In  addition,  the 
multiple  weapon  delivery  under  Independent  passes  remains  to  be  solved. 
Questions  such  as  the  effect  of  the  number  of  weapons  delivered  per  pass  o«*  the 
number  of  passes  on  the  variability  in  Pk  have  not  been  answered.  Finally,  it 
is  conceivable  that  the  general  problem  of  weapon  delivery  be  approached  as  a 
A  single  mathematical  model  In  which  the  distribution  of  the  delivery  error  takes 

a  general  form.  The  different  situations  studied  so  far  as  separate  problems 
'  nay  then  be  recovered  as  Special  cases  to  this  general  problem. 
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